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ABSTRACT 

We analyze the spectral asymmetry of Stokes V (circularly polarized) profiles of an individual network patch 
in the quiet Sun observed by Sunrise/IMaX. At a spatial resolution of 0.15" - 0.18", the network elements 
contain substructure which is revealed by the spatial distribution of Stokes V asymmetries. The area asymmetry 
between the red and blue lobes of Stokes V increases from nearly zero at the core of the structure to values 
close to unity at its edges (one-lobed profiles). Such a distribution of the area asymmetry is consistent with 
magnetic fields expanding with height, i.e., an expanding magnetic canopy (which is required to fulfill pressure 
balance and flux conservation in the solar atmosphere). Inversion of the Stokes / and V profiles of the patch 
confirms this picture, revealing a decreasing field strength and increasing height of the canopy base from the 
core to the periphery of the network patch. However, the non-roundish shape of the structure and the presence 
of negative area and amplitude asymmetries reveal that the scenario is more complex than a canonical flux tube 
expanding with height surrounded by downflows. 

Subject headings: Sun: surface magnetism — Sun: atmosphere — Polarization 



1. INTRODUCTION 

Asymmetries in spectral profiles of Stokes V are 
commonly measured on the solar surface, from the 
penumbrae of sunspots to the quiet Sun; they have 
been known for decades (e.g. Solanki & Ste nflol 1 19841; 
Grossmann-Doerth & Schiissler 1996; Martinez Pillet et al. 
19971; ISigwarthl 120001 1200 U LRezaei et all I20CV 



Martin ez Gonzalez et alj|2008t IViticchie & Sanchez Almeid; . 
20111) . The generally accepted model to produce asymmetric 



Stokes V profiles is by having anticorrelated v elocity and 
magn e tic field gradients along the line of sight (|Illing et aT] 
\wm ISanchez Almeida et al] [T989t lLopez Aristd 120021) . 
Therefore, circular polarization asymmetries encode the 
structuring of magnetic fields along the line of sight. How- 
ever, the limited spatial resolution of most observations mixes 
the contributions of different regions to the asymmetries, 



hence, allowing only indirect inferences of gradients along 
the line of sight. 

Outside of sunspots the classical model to produce asym- 
metries is at the canopies of the magn etic fe aures exp anding 
with height (iGrossmann-Doerth et all 11988b ISolankil fl989h . 
This model has been able to reproduce the simultaneous 
presence of Stokes V are a asymmetry, the absence of a Stokes 
V ze ro-crossing shift (Solanki 1986; Martinez Pillet et aT] 
119971) . and the c hange in sign of t he area asymmetry near 
the solar limb (Biinteetal. 1993) in obervations of the 
quiet Sun a nd active region plage taken at rather low spatia l 
resolution dStenflo et alj 119871; Martinez Pillet et all 11997b . 
Alternatives have also been suggested, such as the pres- 
ence o f bulk veloc i ties in the flux tube dBellot Rubio et all 
Il997allbl : ISigwarthl l200lh. or the MISMA and M ESMA 
models (ISanchez A lmeida & Landi Pegl' Innocent.il 119961 : 
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ICarroll & Kopl2007h . 

This letter is devoted to employing Stokes V asymmetry 
in order to diagnose the internal structure of network mag- 
netic features, i.e., the stongest flux a ccumulations in the quiet 
Sun. It is known, since the work of Stenflo ( 1973), that net- 
work patches are composed of magnetic elements with kG 
fields. Later, two very different models have been proposed 
to explain the observed asymmetries. First, the thin flux tube 
model that assumed a rather vertical magnetic tube expanding 
with height. Second, the MISMA hypothesis that assumes 
an infinite number of magnetic atmospheres in the resolution 
element. In both cases, the cause for the asymmetries is the 
correlation between the distribution of magnetic fields and ve- 
locities. However, the MISMA model is a one-dimensional 
model, having no predictive power for the spatial variation of 
asymmetries in magnetic structures. In order to understand 
how magnetic fields are really organized and to discern be- 
tween these two models (or create new ones), we must ob- 
serve at higher spatial resolutions, which will allow us to di- 
rectly map the substructure of magnetic features and to trace 
gradients along the optical depth. 

2. DATA AND ANALYSIS 

We map the Stokes V profil e asymmetries in the quie t Sun 
using da ta obtained by IMaX dMartmez Pillet et alJl20Tlb on 
Sunrise dSolankiet all 12010: Bart holetalJl2011l) . The data 
analyzed here consist of Stokes / and V profiles of the Fe i 
5250.2 A line (Lande factor of 3) sampled at 1 1 wavelegth 
positions that were recorded on June 10th at the disk center. 
The effective field of view (FOV) is 46.8"x46.8" and the spa- 
tial resolution is about 0.2". The noise level in Stokes V is 
10~ 3 in units of the continuum intensity, I c , prior to the appli- 
cation of the phase diversity reconstruction technique. After 
this procedure, the noise of the reconstructed data is spatially 
and temporaly correlated, but still follows a Gaussian distri- 
bution, with a standard deviation of 2.5 x 10~ 3 I c . 

For the analysis presented in this letter, we select a small 
region of the FOV containing the largest and the most intense 
circular polarization patch (see Fig. [TJ. It has a magnetic flux 
of ~ 2 x 10 19 Mx contained in an area of ~ 5.5 Mm . The pat- 
tern of the granulation is distorted within the magnetic patch, 
which is typical of kG field accumulations. In order to min- 
imize the effect of the noise in our analysis, we select only 
those profiles with amplitudes larger than 5 times the noise 
level. 46.6% of the non-reconstructed, and 29.3% of the re- 
constructed profiles belonging to the small selected area meet 
this criterion. The selected Stokes V profiles in this region are 
mostly regular profiles. More precisely, 93.6% and 88.9% of 
the profiles have two lobes in the non-reconstructed and in the 
phase diversity reconstructed data, respectively. We define the 
area (SA) and amplitude (6a) asymmetries of regular profiles 
as: 
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The symbols A t = 5250.0164 A and A f = 5250.4014 A are 
the initial and final wavelengths of the integral over the ob- 
served profile. The symbol <x is the sign of the Stokes V blue 
lobe. The ab and a r symbols are the unsigned blue and red 
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Figure 1. From top to bottom, continuum intensity, circular polarization 
amplitude and full width at half maximum of the intensity profile for the 
stongest flux region in the IMaX field of view. Left panels display the ob- 
servables of the original data, while the right ones show the same scenes after 
phase-diversity reconstruction. The white contours represent the isocontours 
of log[Av// c ] = -2.3,-1.9 for the non recostructed and the reconstructed 
data, respectively. The symbol Ay represents the amplitude of the Stokes V 
profile. 

lobe amplitudes, respectively. The non-linearity of the equa- 
tions to compute asymmetries make it difficult to evaluate the 
relevance of instrument errors. The first error source to take 
into account is the photon noise. We have verified (simulating 
Stokes V profiles as derivatives of a gaussian that match IMaX 
intensity data and including gaussian noise) that the photon 
noise does not introduce a bias in the calculation, only spreads 
the results, the relative errors in the area and amplitude asym- 
metries being at most ±15%. 

Calibrations performed at the launch station in Kiruna, 
proved that there exist considerable amounts of crosstalk be- 
tween linear and circu lar polarization introduced by the ISLiD 
dGandorfer et al.l201 ll) . The relative amounts of Stokes Q and 
Stokes U change with time due to image rotation and cannot 
be predicted for a specific data set. When IMaX observes 
in a vectorial mode, all four Stokes parameters are measured 
and multiplication by the inverse of the Mueller matrix takes 
into account all crosstalk sources (leaving as unknown only 
the zero of the azimuth angle). I n the longitudinal observing 
mode used here (L12-2, refer to Martinez Pillet et all 1201 lb 
only Stokes / and Stokes V are measured making a proper cor- 
rection of crosstalk impossible. However, we expect the im- 
pact of this crosstalk to be negligible for the investigation per- 
formed in this work that concentrates on a relatively large net- 
work patch. First of all, the fields in these regions are known 
to be predominantly vertical and all horizontal fields observed 
with IMaX in other full vector maps show that the linear polar- 
ization signals mainly corresponded to internetwork elements 
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dDanilovic et al.l2010l) . Second, inspection of similar network 
patches obser ved in the IMaX vect or mode V5-6 (bottom part 
of Figure 3 in lSolanki et al. 2010), show that they never dis- 
play linear polarization signals above 3 sigma levels at the 
network concentration, neither inside it nor at the boundaries, 
where a possible canopy could generate some linear polariza- 
tion. We thus conclude that in the region under study, the 
known crosstalk in the longitudinal mode of IMaX can be ig- 
nored even though it cannot be accounted for as is done in the 
IMaX vector modes. 

Finally, the effect of the finite spectral (and spatial) resolu- 
tion, as well as the effect of the phase-dive rsity reconstruction, 
is studied in lAsensio Ramos et al.l d201 lb for the special case 
of MHD simulations. These authors conclude that the spatial 
degradation tends to understimate both amplitude and area 
asymmetries. Degrading the spectral resolution, however, 
tends to oveerstimate the area asymmetry while it does not 
affect str ongly the amplitude asym metries, confirming earlier 
results of Solanki & Stenflo ( 1986). After applying the recon- 
struction algorithm, the computed asymmetries resemble the 
original ones, i.e., for MHD simulations at full spectral and 
spatial resolution. However, we must note that these conclu- 
si ons depend on the actual shape of the profiles. The work 
of Solanki ( 1986) predicts (for the IMaX case) an overstimate 
of the area asymmetries by a factor ~ 1.6 and an understi- 
mate of the amplitude asymmetries by a factor ~ 0.8. Conse- 
quently, strongly asymmetric/antisymmetric profiles observed 
by IMaX do correspond to profiles with large/small asymme- 
tries at full spectral resolution, although the exact values may 
differ. 

3. CIRCULAR POLARIZATION ASYMMETRIES AND THE 
MAGNETIC FIELD GEOMETRY 

Figure |2] shows maps of the area and amplitude asymme- 
tries of the observed Stokes V profiles in the selected patch. 
The left panels correspond to the non-reconstructed data and 
the right ones to the phase-diversity reconstructed ones. All 
four panels display the same pattern: the smaller asymmetries 
tend to be located in the inner parts of the structure, while the 
larger ones lie at the borders, close to where the signal drops 
below the noise level. The larger asymmetries situated at the 
borders of the structure have dominantly positive values, but 
negative values are also found (e.g. along the top edge of the 
feature). The spatial pattern for the two types of asymmetries 
is very similar, the correlation coefficient between amplitude 
and area asymmetries being 0.9. 

Figure [3] displays three sets of Stokes / and V profiles, ex- 
tracted from the positions marked with small grey arrows in 
Fig. [2] The Stokes V profile displayed in the left panels 
(6 A - 0.5% and 6a - 4%) is representative of the nearly 
antisymmetric profiles found in the central part of the patch. 
The profile in the middle panels has larger area and amplitude 
asymmetries (38% and 33%, respectively), as found closer 
to the boundary of the magnetic feature. Finally, the pro- 
file in the rightmost panels is a one-lobed profile, with ex- 
treme asymmetries: 81% area asymmetry, and 80% amplitude 
asymmetry. Such profiles are found almost exclusively at the 
very edge of the feature. 

In order to check if the expanding flux-tube model applies 
to our observations, we use the SIRJUMP cod e, which is a 
modifi cation of the SIRGAUS code described in lBellot Rubiol 
( 20031). The two codes follow the inversion strategy of SIR 
( Ruiz Cobo & del Toro Iniestall 19921) . SIRJUMP has recently 
been employed to interpret single-lobed Stokes V profiles in 



the quiet Sun bv lSainz Dalda etal.1 (120121) . With SIRJUMP 
we retrieve the parameters of a model atmosphere showing 
discontinuous stratifications of some of the parameters along 
the line of sight. The position of the jump (i.e., the discontinu- 
ity) is a free parameter in the inversion scheme. The code can 
put it below optical depth unity, which would indicate that 
no discontinuity is needed to explain the observed profiles. 
Above and below the jump, the magnetic field and the line-of- 
sight velocity are assumed to be constant with optical depth. 
The temperature stratification is modified with respect to an 
initial guess by the inversion algorithm in three nodes and in- 
terpolated between them. The microturbulent velocity and the 
inclination of the magnetic field vector to the line of sight are 
height independent. Since we do not have the information of 
the linear polarization, we cannot retrieve the magnetic field 
azimuth, and can barely infer the inclination. Therefore, we 
perform two different inversions, one forcing the field to be 
aligned with the line-of-sight and another one leaving the in- 
clination as a free parameter. The macroturbulent velocity is 
fixed to 2.06 km s _1 , which matches the broadening induced 
by the Point Spread Function of the IMaX instrument. 

Figure |4] displays the results of the two different inversions. 
Top left panels show the optical depth at which the discon- 
tinuity occurs for the inversion with the inclination as a free 
parameter (panel 4a) and the inversion forcing the inclination 
to 0° (panel 4b). In both cases, the results are very similar, 
i.e., they do not depend much on the assumed geometry of 
the magnetic field. For logTsooo ^ 0.5 the contribution to the 
spectral line is small and the discontinuity is not well con- 
strained. In spite of the noise some trends can be seen. In the 
central part of the structure we retrieve a magnetic field with 
no clear discontinuity. As we move from the center towards 
the edges of the structure, the discontinuity appears at in- 
creasingly higher layers. Therefore, the inversions made pixel 
by pixel have naturally retrieved the expansion of a magnetic 
flux concetration with height in the solar atmosphere. At all 
heights (below and above the discontinuity), the pattern of the 
velocity is qualitatively similar, following the local distorted 
granulation. 

In the lower panels of the same figure we display the rel- 
ative differences of the longitudinal components of the mag- 
netic field (panels e and f) and the velocity (panels g and h) 
above and below the discontinuity. Note that, since we do not 
have linear polarization, the only trustable quantity is the lon- 
gitudinal component of the magnetic field. We have used the 
following definitions: 
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where the superindex a and b refer to the quantities above and 
below the discontinuity, respectively. From these figures, we 
can see that both the longitudinal magnetic field and the line- 
of-sight velocity are almost constant with optical depth at the 
core of the structure, coinciding with strong kG fields. As we 
move towards the borders of the structure (where the longitu- 
dinal component of the magnetic field is of the order of 100 
G), the relative variation of both quantities increases. Note 
that the gradient of the longitudinal magnetic field changes 
sign at the regions where the area and amplitude asymmetries 
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Figure 2. Top panels represent the area asymmetry of the circular polarization profiles of the selected strong flux patch. Bottom panels display the amplitude 
asymmetry. Right (left) panels present the results showing the reconstructed (non-reconstructed) data. The areas in white represent those pixels where the signal 
to noise ratio of the Stokes V profiles is lower than 5. Grey arrows are the locations of the profiles displayed in Fig. [3] 

pattern, in the IMaX FOV. The core of the structure is, on 
average, characterized by an area asymmetry of 0.2% and an 
amplitude asymmetry of 4.4% (with an upper limit for the er- 
ror of +1.3% for both averaged area and amplitude asymme- 
tries). In the outermost parts of the network patch, the Stokes 
V profiles have the largest area and amplitude asymmetries, 
even aproaching 100% (one-lobed profiles). The fact that we 
are seeing a spatial variation of Stokes V asymmetries means 
that, for the first time, we are resolving the internal structure 
of a magnetic patch belonging to the photospheric network. 
The non negligible amplitude asymmetry inside the magnetic 
feature suggests that unresolved motions are still present (def- 
initely gradients along the line of sight and probably horizon- 
tal gradients). 

We inverted our data using the SIRJUMP code, which mod- 
els the gradients of magnetic field and line-of-sight velocity 
with an abrupt discontinuity at a certain optical depth. From 
two different assumptions on the inclination of the field, we 
inferred that the observed network patch is made of resolved 
kG fie lds which is consistent with the work of lLagg et alJ 
(2010). Moreover, the inversions naturally retrieve a geom- 
etry corresponding to a magnetic field expanding with height. 

The spatial distribution of the Stokes V asymmetries and 
the inversions applied to these data suggest a geometry simi- 
lar to t he classical model of magne tic features expanding with 
height (Sp ruitll976tlSolankill993h that takes pressure balance 
and magnetic flux conservation with height into account and is 
consistent w ith three-dimensional magn eto-hydrodynamical 
simulations (I Ye lies Chaouche et al.l 120091) . In such a model, 
the central core of the structure traces the body of the flux 
tube, containing the vertical kG fields, while the external parts 
correspond to the canopy, i.e. expanding, inclined field over- 
lying the quiet photosphere (this expansion is caused by the 
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Figure 3. Examples of Stokes V profiles. The left panels display a respre- 
sentative of the relatively antisymmetric profiles found in the center of the 
structure, middle panels the more asymmetric ones found near the borders 
of the patch, while the right panels depict a one-lobed profile at the edge of 
the magnetic patch. The exact locations are marked in Fig. |2]by small grey 
arrows. The observed profiles are represented by the rhombs and the best fit 
using an inversion algorithm based on a discontinous model along the optical 
depth (see the text) is displayed with a blue line. The inversion with the incli- 
nation forced to zero is represented by the green lines and the inversion with 
the inclination as a free parameter is displayed with the blue lines. 



were negative. 



4. CONCLUSIONS 



We have computed the area and amplitude asymmetries of 
the Stokes V profiles of the Sunrise/IMaX dataset with the 
best spectral sampling (35 mA). We selected the strongest 
magnetic flux accumulation, which belongs to the network 
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Figure 4. Results of the two inversions of the original data based on a model that has a discontinuity along the optical axis. Panels a, c, e, and g display the 
inferred quantities for the inversion in which the inclination of the magnetic field is a free parameter. Panels b, d, f, and h represent the retrieved parameters from 
the inversion in which we have fixed the inclination of the magnetic field with respect to the line-of-sight to 0° . Similar results can be derived from the phase 
diversity reconstructed data. Panels a-d show the optical depth of the discontinuity (left) and the temperature at logrsooo = — 1.5 (right). Panels e-h display the 
relative differences of the line-of-sight magnetic field (left) and velocity (right) above and below the discontinuity (see Eq. 3) 



need to conserve magnetic flux for decreasing gas pressure 
with height). 

The standard model for the creation of asymmetries in the 
presence of such a canopy is based on the presence of a dis- 
continuity in both the magnetic field strength and the line- 
of-sight velocity at the base of the canopy (i.e. abrupt gra- 
dients in such quantities along the optical depth). In the 
past, it was assumed that there is no velocity flow above the 
canopy's base but there is a downflow below it, with the mag- 
netic field being present only above the canopy base. This ge- 
ometry produced a positive area (and amplitud e) asymmetry 
dGrossmann-Doerth et all 119881: ISolankil 119891) . wit h nearly 
one-l obed profiles in the outer parts of the flux tubes (ISolankil 
1989) in agreement with the observations. 

Our results for the asymmetries (see Fig. |2} are in agree- 
ment with this general picture. However, there are some lo- 
cations along the edge of our feature where the area (and 
amplitude) asymmetry is negative, indicating that either the 
magnetic field or the velocity field gradient must change sign 
there. We know that the kG network fields are embedded in 
the more disorganised weakly magnetized quiet photosphere. 
Hence, we would expect a smooth transition between the or- 
ganized kG fields to the tangled weak photospheric fields 
(note that we have selected the profiles above the noise level, 
i.e., the borders of the structure depend on the amplitude rec- 
quired for the profiles). If those regions are mostly tracing 
the quiet photosphere (note that the height of the discontinu- 
ity is located at lower layers and that the temperature is colder 
than that of the central part of the structure), it is not sur- 
prising to find gradients of magnetic field and velocity that 
change sign as compared to the classical canopy overlying a 
non-magnetized photosphere. 

The observations presented in this letter support the inter- 
pretation of the network as kG fields expanding with height. 
In contrast to what was assumed in the past, we find veloc- 
ity patterns above and below the canopy that follow the local 
distorted granulation. To strengthen the model for strong flux 



concentrations as the one studied in this letter, in the future 
we need to include in our analysis the constraints introduced 
by the linear polarization and, more importantly, the simulta- 
neous information of many spectral lines to trace the canopy 
at different heights. 
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